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We report on a search for anomalous production of Z boson pairs through a massive resonance 
decay in data corresponding to 2.5-2.9 fb~^ of integrated luminosity in pp collisions at -^s = 
1.96 TeV using the CDF II detector at the Fermilab Tevatron. This analysis, with more data 
and channels where the Z bosons decay to muons or jets, supersedes the 1.1 fb~^ four-electron 
channel result previously published by CDF. In order to maintain high efficiency for muons, we use 
a new forward tracking algorithm and muon identification requirements optimized for these high 
signal-to-background channels. Predicting the dominant backgrounds in each channel entirely from 
sideband data samples, we observe four-body invariant mass spectra above 300 GeV/c^ that are 
consistent with background. We set limits using the acceptance for a massive graviton resonance 
that are 7-20 times stronger than the previously published direct limits on resonant ZZ diboson 
production. 



PACS numbers; 12.60.Cn,13.85.Rm,14.70.Hp 

I. INTRODUCTION 

The standard model of particle physics (SM) 
has been enormously successful, but many key 
questions remain to be answered by a more 
complete theory. New theoretical ideas can be 
tested with collider experiments, but it is also 
worthwhile for experiments to search broadly for 
anomalous "signatures" . One common tactic 
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is to look where experiments are keenly sensi- 
tive. The consummate example of this method is 
the Z' boson search, in which a low-background, 
well-understood observable (the dilepton invari- 
ant mass) is used to constrain the new physics 
models that predict a dilepton resonance. Dibo- 
son resonance searches are an attractive analog 
of the Z' boson searches, involving higher mul- 
tiplicities of the same outgoing particles and ad- 
ditional mass constraints, both of which serve to 
further suppress experimental backgrounds. Di- 
bosons are the dominant channels for high mass 
higgs searches, and new physics scenarios pre- 
dict particles such as Randall-Sundrum gravitons 
which would decay into dibosons [1]. The irre- 
ducible SM diboson background processes occur 
at such a low rate that they have only recently 
been observed at the Tevatron [2, 3] at low dibo- 
son mass (Mzz < 300 GeV/c^). At high diboson 
mass {Mzz > 300 GeV/c^) the total backgrounds 
are tiny. 

This article presents a search for a diboson res- 
onance in data corresponding to 2.5-2.9 fb~^ of 
integrated luminosity in ^/s = 1.96 TeV pp col- 
lisions at the CDF II detector at the Fermilab 
Tevatron, in the decay channel X — >■ ZZ. These 
ZZ diboson processes have been well-studied at 
the LEP experiments, which observed no signif- 
icant deviation from the standard model expec- 
tation up to an e^e~ center-of-mass energy of 
207 GeV [4-7]. The LEP data place only in- 
direct constraints on heavier, resonant ZZ di- 
boson production [8], however, and direct pro- 
duction constraints at high ZZ diboson masses 
must be probed at hadron colliders. To the pre- 
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viously published CDF search for four-electron 
production via X — > ZZ — >■ eeee with data cor- 
responding to 1.1 fb^^ of integrated luminosity 
[9], we now add the dijet channels eejj and ^i^ijj, 
which improve sensitivity at very high X masses 
where their background is negligible, and the 
four-electron or -muon channels ee/i/i and (XfJ-fi^i, 
which contribute sensitivity to new physics at in- 
termediate masses where the Z + jet(s) back- 
grounds are larger. 

Because there are four or more outgoing lep- 
tons or quarks, the analysis is sensitive to sin- 
gle lepton and jet reconstruction efficiencies to 
approximately the fourth power. In particular, 
events may often have one or more leptons with 
1 7^1 > 1 where muon acceptance and tracking ef- 
ficiency are lower than for \r]\ < 1. Consequently 
the four-lepton channels motivate development 
and use of techniques to improve electron and 
muon reconstruction and identification efficien- 
cies while exploiting the kinematics of the sig- 
nature to keep backgrounds low. To augment 
forward muon coverage, the present analysis also 
employs, for the first time, a new method of re- 
constructing charged particles in the silicon de- 
tectors using constraints from particle traces in 
forward regions with partial wire tracker cover- 
age. 

The aim of the search is sensitivity to any mas- 
sive particle that could decay to ZZ. Though we 
avoid focus on any one specific model, we choose 
a benchmark process that is implemented in sev- 
eral popular Monte Carlo generation programs, 
the virtual production of gravitons in a simple 
Randall-Sundrum RSI scenario [1, 10], to fix ac- 
ceptance for the search and quantify its sensi- 
tivity. The geometry of the model consists of 
two three-branes separated from each other by 
a single extra dimension. Boundary conditions 
at the branes quantize the momentum in the ex- 
tra dimension, leading to a Kaluza-Klein tower 
of discrete, massive gravitons. In RSI scenarios 
with the standard model particles confined to ei- 
ther brane, a discovery would involve graviton de- 
cays directly to photons or leptons — the graviton 
branching ratio to ZZ is significant but the Z bo- 
son branching ratio to leptons is small. In more 
complex but well-motivated scenarios with SM 
particles allowed in the extra dimension, how- 
ever, graviton decays to photons, leptons, and 
light jets can be suppressed, and dibosons be- 
come an important discovery channel [11, 12]. 

The organization of this article is as follows: 
Section 2 describes relevant components of the 



CDF II detector; Section 3, the event selection; 
Section 4, the background estimates; and Section 
5, the results. 



II. THE CDF II DETECTOR 

The CDF II detector is a general purpose mag- 
netic spectrometer surrounded by electromag- 
netic and hadronic calorimeters and muon de- 
tectors designed to record Tevatron pp collisions. 
We briefly describe the components of the detec- 
tor relevant to this search. A complete descrip- 
tion can be found elsewhere [13]. 

A combination of tracking systems recon- 
structs the trajectories and measures momenta of 
charged particles in a 1.4 T solenoidal magnetic 
field. Trajectories of charged particles are recon- 
structed using an eight-layer silicon microstrip 
vertex tracker [14] at radii 1.3 < r < 29 cm from 
the nominal bcamlinc^ and a 96-layer open-cell 
drift chamber (COT) providing eight superlay- 
ers of alternating stereo and axial position mea- 
surements [15] at large radii 43 < r < 132 cm. 
The COT provides full geometric coverage for 
\r]\ < 1.0. The average radius of the lowest ra- 
dius axial (stereo) COT superlayer is 58 (46) cm, 
providing partial coverage for \r]\ < 1.7(1.9). The 
silicon tracker provides full coverage for jyy] < 1.8. 

Outside the tracking volume, segmented elec- 
tromagnetic (EM) lead-scintillator and hadronic 
(HAD) iron-scintillator sampling calorimeters 
measure particle energies [16]. The central {\r]\ < 
1.1) calorimeters are arranged around the in- 
teraction point in a projective-tower cylindrical 
geometry, divided azimuthally into 15° wedges. 
This calorimeter measures electron energies with 
a resolution of [a{E)/E]^ = {13.5%)^ /Et + 
(2%)^. The forward calorimeters (1.1 < \r]\ < 
3.6) are arranged in an azimuthally-symmetric 
disk geometry and measure electron energies with 
a resolution of [a{E)/E]^ = (16.0%) V^J-H (1%)^. 
Wire chambers (scintillator strips) embedded in 
the central (forward) EM calorimeters at ~ GXq, 



^ CDF uses a cylindrical coordinate system in which 9 
{<j>) is the polar (azimuthal) angle, r is the radius from 
the nominal beam axis, and +z points along the proton 
beam direction and is zero at the center of the detec- 
tor. The pseudorapidity is defined as rj = — lntan{0/2). 
Energy (momentum) transverse to the beam is defined 
E'Y = -Bsinfi (px = psinS), where E is energy and p is 
momentum. 



6 



the average depth of shower maximum, provide 
position and lateral shower development mea- 
surements for 1 77 1 < 2.5. 

Beyond the calorimeters, muon drift chambers 
and scintillators measure particles that traverse 
the entire inner and outer detectors and reject the 
instrumental backgrounds of the central muon 
triggers. The central muon chambers (CMU) 
lie just outside the central hadronic calorimeter 
with (^-dependent coverage for 0.03 < \ri\ < 0.63. 
The central muon upgrade (CMP) augments the 
CMU coverage in (p and lies behind another ap- 
proximately 3 interaction lengths of steel. The 
central muon extension (CMX) extends coverage 
into the region 0.65 < < 1.0. 

The beam luminosity is determined by mea- 
suring the inelastic pp collision rate with gas 
Cherenkov detectors [17], located in the region 
3.7 < |?7| < 4.7. 

At each bunch crossing, a three-level trig- 
ger system [13] scans the detector output for 
\r]\ < 1.1 electrons or I77I < 1.0 muons with at 
least 18 GeV/c of transverse momentum. We ac- 
cept events that satisfy one of four trigger paths: 
one that requires a deposition of at least 18 
GeV transverse energy in the calorimeter con- 
sistent with an electron and a matching COT 
track with at least 9 GeV/c of transverse mo- 
mentum; another with fewer electron identifica- 
tion requirements intended to ensure high effi- 
ciency for very energetic electrons; a muon path 
requiring a COT track with at least 18 GeV/c 
of transverse momentum pointing toward signals 
in both the CMU and CMP chambers (a CMUP 
trigger) and traversing the calorimeter consistent 
with a minimum-ionizing particle; or a similar 
muon path with signals in the CMX chamber in- 
stead of the CMU and CMP chambers. 



III. DATA COLLECTION 

We use data corresponding to an integrated lu- 
minosity of 2.5-2.9 fb~^ depending on the data 
quality criteria applicable to the relevant ZZ di- 
boson decay channel. We separately analyze six 
channels: eeee, eeyit/Lt, niiee, /x/z/z/x, eejj, and 
^ij-iii- Events are divided into the six categories 
based on the trigger, where the first lepton de- 
notes the required trigger path, and the presence 
of lepton and jet candidates identified using the 
criteria listed in Tables I through IV. The trig- 
ger lepton criteria are the most stringent; sub- 
sequent kinematic signature selections yield very 



TABLE I: Calorimeter electron identification criteria. 
We require Had/EM, the ratio of energies measured 
in the hadronic and electromagnetic calorimeters, to 
be less than f{E) = 0.055 + 0.00045 x [E/GS^) 
where E is the measured calorimeter energy. Isocai 
is the calorimeter energy measured within Ai? = 0.4 
centered on the electron, excluding the electron en- 
ergy, r^dct is calculated assuming an origin at 2 = 0. 
LshrTrk is a lateral shower shape variable described 
in Ref. [18]. 



Criteria 


Trigger 


Central 


Forward 


£;T(GeV) 


> 20 


> 5 


> 5 


Track zq (cm) 


< 60 


< 60 




Had/EM 


<f{E) 


<f{E) 


< 0.05 


ISOcal/ Ecal 


< 0.2 


< 0.2 


< 0.2 


\mct\ 






< 2.5 


LshrTrk 


< 0.4 






Track pT (GeV/c) 


> 10 







low backgrounds, allowing very efficient identi- 
fication criteria to be used for the other lepton 
candidates. Events accepted by either electron 
trigger path and containing at least one electron 
candidate that fired the electron trigger and sat- 
isfied the offiine selection criteria are excluded 
from the muon-triggered categories. There are 
no events that satisfy the requirements of more 
than one category. 

During this selection we identify the events 
containing at least two leptons (including the 
trigger lepton) using the nominal CDF event re- 
construction software. We reprocess these events 
using a revision of the software with improved 
tracking, including more efficient forward track- 
ing algorithms, and then select all final particles 
from the reprocessed data. In this way, we avoid 
CPU-intensive reconstruction of a 96.6% subset 
of the sample that has no chance to pass our fi- 
nal selection. Nevertheless, two subsets of the 
data corresponding to 200 pb^^ of integrated lu- 
minosity each were fully reprocessed without the 
initial two-lepton selection and analyzed to con- 
firm that the procedure used for the full dataset 
is fully efficient for events of interest to the anal- 
ysis. 

The electron criteria listed in Tables I through 
II are nearly identical to the previous eeee analy- 
sis [9] . The double Z boson mass peak signature 
admits little background, allowing appreciably 
more efficient electron criteria than those used 
for many other CDF analyses (for example, Ref. 
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TABLE II: Track electron identification criteria. 
Tracks must consist of measurements in several COT 
superlayers. Silicon measurements are not required. 
Isotrk is the scalar sum of the momenta of all tracks 
measured within a circle of AR — 0.4 centered on the 
electron track direction. ARem is the separation in 
the rj — (f) plane between the electron track and the 
nearest calorimeter electron cluster, as defined in Ta- 
ble I. 

Criteria 

PT (GeV/c) > 10 

Axial Superlayers > 3 

Stereo Superlayers > 2 

I Track zo\ (cm) < 60 

Ptrk / (ISOtrk + Ptrk) > 0.9 

, , , I 200 um silicon 

2 mm no silicon 
ARem > 0.2 



[13]). We require either an isolated calorimeter 
cluster with electron-like energy deposition or, to 
recover acceptance, an isolated track pointing at 
uninstrumented regions of the calorimeters. Such 
regions constitute approximately 17% of the solid 
angle for {rjl < 1.2 and would otherwise reduce 
our four-electron acceptance by a factor of two. 
The transverse energy threshold for non-trigger 
electrons is 5 GeV. As the mass of the signal res- 
onance X — >■ ZZ increases, the energies of the 
two Z boson decay products become asymmetric 
in the detector frame, and thus our criteria must 
efhciently select leptons with transverse momen- 
tum of order 10 GeV/c as well as leptons with pT 
of hundreds of GeV/ c. 

The muon criteria listed in Table III require 
an isolated track satisfying basic track quality 
criteria and depositing minimal energy in the 
calorimeter. We make use of a new track recon- 
struction algorithm, described in the Appendix, 
and apply less stringent energy and isolation re- 
quirements than typical for CDF high pT analy- 
ses so as to increase our acceptance and efficiency. 
Muon candidate tracks may be matched to in- 
formation in the muon chambers, but to recover 
acceptance lost due to gaps in chamber cover- 
age and to recover efficiency lost due to pointing 
requirements, chamber matching is not required 
except for the trigger muon. 

Jets must satisfy the criteria listed in Table IV 
and jet energies are corrected for instrumental ef- 



TABLE III: Muon identification criteria. The CMU, 
CMP, and CMX match variables compare the track 
position extrapolated to the relevant muon chambers 
with the chamber position measurements. The non- 
trigger pt requirement is lower for tracks with muon 
chamber information attached. IsOcai is the sum of 
calorimeter energies measured in towers within a cir- 
cle of A_R = 0.4 centered on the muon tower. Eem 
and Ehad are the electromagnetic and hadronic 
calorimeter energies recorded in towers intersected by 
the muon track, and /em(p*'^'^) = 4 + max(0, 0.0115* 
((fsVT^-lOO)) and/^fAub*"'') = 12 + max(0, 0.028* 

((fevTc ~ lOO)) are functions of the track momentum. 
The cuts on track curvature k, its uncertainty o^., 
and the probability of the fit Prob(x'^, Jidof) reject 
poorly measured tracks. 



Criteria 


CMUP 


CMX 


Non-trigger 


p^''=(GeV/c) 


> 20 


> 20 


> 2,10 


CMU match 


< 10 cm 






CMP match 


< 20 cm 






CMX match 




< 10 cm 





common to all categories 



Iso'=^7p*'^'= < 0.2 

EEM{Ge^) < /iJA/(p*''') 



EHADiGeV) 


< fHADip''-'') 




> 2.5 


Prob(x2,7idof) 


> 10-10 


\zo\{cm) 


< 60 


\do\ 


J 200 /im silicon 


1 2 mm no silicon 



fects [20]. Before relying on the jet energy mea- 
surements for the two-lepton two-jet analysis, we 
have verified that these corrections balance trans- 
verse momentum in the Z-l-jet(s) events consid- 
ered here. We choose kinematic requirements on 
individual jet energies, on dijet invariant masses, 
and on four-body masses involving jets so that 
the systematic uncertainties on X ^ ZZ signal 
acceptances and efficiencies from mis-modeling of 
QCD radiation or jet reconstruction effects are 
small. 

Figs. 1 and 2 show comparisons of the peak- 
ing and background components of the dielec- 
tron and dimuon yield. The combination of our 
changes to the identification criteria and to the 
tracking algorithms increases the peak yield by 
factors of 1.8 and 4.3, respectively. 

Measurements of the pp 'y*/Z — > ee and 
pp — > ^* /Z — > /i/i cross-sections provide an im- 
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TABLE IV: Jet identification criteria. The JETCLU 
algorithm is discussed in Ref. [19]. ARem is the 
separation in the rj — <j) plane between centroids of 
the jet cluster and the nearest electron cluster. 

Selection Criteria 



Algorithm 
£;5?^(GeV) 

I ^centroid | 
EEM/Etot 

ARem 



JETCLU 0.4 Cone 

> 10 

< 3.64 

< 0.95 

> 0.4 




£ = 2.9 fb' 
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FIG. 2: Z — >■ fifi yield and background comparison 
between dimuon candidates consisting of a trigger 
muon and an muon candidate satisfying either (lower 
set of points) the CDF standard muon criteria or (up- 
per set of points) the criteria employed in the present 
analysis. The peak yield increases from about 35,000 
candidates with the standard criteria and tracking to 
150,000 candidates with our optimized criteria and 
tracking. The continuum background increases by 
a factor of 14 for 81-101 GeV/c^. This background 
is later suppressed by the four-body kinematic selec- 
tion. 



FIG. 1: Z — > ee yield and background comparison 
between dielectron candidates consisting of a trigger 
electron and an electron candidate satisfying either 
(lower set of points) the CDF standard electron cri- 
teria or (upper set of points) the criteria employed in 
the present analysis. The peak yield increases from 
about 146,000 with the standard criteria to 256,000 
candidates with our optimized criteria. The corre- 
sponding increase in continuum background is mod- 
est, a factor of 2.0 for 81-101 GeV/c^. This back- 
ground is later suppressed by the four-body kine- 
matic selection. 
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portant test of our understanding of the trigger, 
reconstruction, and identification efficiencies and 
Monte Carlo modeling for this new lepton selec- 
tion. We divide the data into 18 data-taking pe- 
riods and measure each of the above efficiencies 
for each period. For each period, we then com- 
pute the Drell-Yan cross-section for all combina- 
tions of trigger and lepton type using a signal 
plus background fit to the data and Drell-Yan 
Monte Carlo. The average instantaneous lumi- 
nosity tends to increase with data-taking period 
as Tevatron upgrades were brought online. Figs. 
3 and 4 show the resultant cross-sections and 
their dependence on time. 

For the cases where both electron energy mea- 
surements come from the calorimeters, the fit sig- 



FIG. 3: Z — > ee cross-sections (and averaged cross- 
sections) for 66 < Mu < 116 GeV/c^ and various 
selections. The horizontal axis indicates the 18 peri- 
ods for five selections in succession: two trigger elec- 
trons (252.1±1.2 pb), a trigger electron and a central 
calorimeter electron (TRIG-(-CEM, 248 ± 1.1 pb), a 
trigger electron and a forward calorimeter electron 
(TRIG-I-PEM, 246.2±0.9 pb), a trigger electron and 
a track electron (TRIG+TRACK, 262.1 ± 2.3 pb), 
and, calculated separately, the combination of a trig- 
ger electron and any electron selected using the anal- 
ysis criteria (249.4 ± 1.6 pb). The averaged cross 
sections are indicated by horizontal lines. Uncertain- 
ties are statistical only with the correlated luminosity 
uncertainty not shown. 
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FIG. 4: Z — )■ /i/i cross-sections (and averaged cross- 
sections) for 66 < Mu < 116 GeV/c^ and various 
selections. The horizontal axis indicates the 18 pe- 
riods for two selections in succession: a CMUP trig- 
ger muon and another muon satisfying the analysis 
criteria (250.0 ± 1.2 pb) and a CMX trigger muon 
combined with another selected muon (263.4 ± 1.8 
pb). The averaged cross sections are indicated by 
horizontal Hnes. Errors are statistical only with the 
correlated luminosity uncertainty not shown. 



nal function is a Breit-Wigner fixed to the world 
average Z boson mass and width convolved with 
a Gaussian resolution function. In the trigger 
plus track lepton sample, where there is a sig- 
nificant radiative tail below the Z boson mass, 
the signal function is a Crystal Ball function. 
In the combined sample both functions are in- 
cluded. In all cases, the background function is 
an exponential. The combined result for the elec- 
tron selection, 249.4 ± 1.6 pb for the mass range 
66 < Mil < 116 GeV/c^, is in agreement with the 
dedicated CDF measurement a^i^/z ^ Br(pp — > 
7*/Z ^ U) = 254.9 ± 3.3(stat.) ± 4.6(sys.) [13] 
and is evidence that the constituent efficiencies 
and scale factors are understood well. 

Note that there are large variations in the trig- 
ger plus track electron cross-section. In these 
cases, the fit to the data underestimates or over- 
estimates the amount of power-law background 
contamination of the low dielectron mass radia- 
tive tail of the Z boson. The independent fit 
to the combination of all selections has a higher 
signal to noise ratio and is not sensitive to this 
effect. 

The results for the muon selection also show 
good agreement with the other CDF measure- 
ments, except for a 20% shift in the first pe- 
riod and up to a 10% shift in many of the later 
periods. The disrepancies are due to imperfect 
modeling of tracking-related muon identification 
efficiencies for the new reconstruction software. 
Based on these results, we assign a 20% system- 
atic uncertainty on the signal acceptance for each 



of the six channels. In part because much of the 
data were collected with the CMUP trigger and 
during periods with measured cross-sections in 
good agreement with expectation, and because 
two of the six analysis channels do not involve 
muons, this systematic uncertainty over-covers 
the observed variation in cross section and is the 
dominant uncertainty for the analysis. Neverthe- 
less, the final sensitivity of the analysis improves 
substantially on the earlier eeee search [9]. 



IV. KINEMATIC ANALYSIS 

After selecting electrons, muons, and jets, we 
consider all possible four-lepton Mi or two- 
lepton two-jet iijj combinations for each event 
that contains a trigger lepton. No requirement 
is made on the mass or charge of dilepton pairs. 
Any two particles must have a minimum separa- 
tion Ai? = \/~Krf^^~K(f? of 0.2. Dilepton pairs 
with tracks present for both Icptons must point 
back to the same zq production location in order 
to suppress background from additional pileup 
interactions. Track timing information is used 
for a very pure veto of muon and track electron 
pairs consistent with cosmic rays. There are no 
events that appear in more than one ZZ diboson 
decay channel. 

For the four-lepton channels, we consider all 
possible combinations of leptons for each event 
and select the one that minimizes a vari- 
able quantifying consistency between the dilep- 
ton masses and the Z boson pole mass: 



Xzz — 



(a4'^ -91.187GeV/c2 



i=l,2 



where (7jv^(i) is the detector mass resolution com- 
puted from the individual lepton calorimeter or 
tracking measurements for the dilepton mass 

and a Gaussian approximation with rrr = 
3.25GeV/c^ allows for the nonzero width of the 
Z boson resonance. 

For the ttjj channels, we consider all possible 
combinations of leptons and jets. We select the 
two highest jets and the dilepton pairing that 
minimizes the first term of the equation above. 
This explicitly avoids possible Z boson mass bias 
in the dijet mass spectrum that would complicate 
the background estimate discussed in the follow- 
ing section. We then require Mz > 20 GcV/c^ for 
each pairing and, for the dijet channels, Xz < ^5 
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for the leptonic Z boson. 

A 'priori we define our signal region to be 
Mx > 300 GeV/c^ so as to avoid most standard 
model backgrounds. For the ££££ modes we fur- 
ther require Xzz ^0 ^^^^ for the ££jj modes we 
require 65 < Mjj < 120GeV/c^. Each event may 
contain additional leptons, jets, or other particles 
beyond the four that contribute to the signal can- 
didate. 
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V. BACKGROUND ESTIMATES 

For both the four-lepton and the dijet chan- 
nels, the dominant backgrounds at high Mx 
are a mixture of Z + jets, + jets, multi- 
jets, and various lower-rate processes resulting 
in one or more hadrons that mimic an electron 
or muon. The diboson processes W^Z — ^ jj££, 
ZZ ££££, and ZZ ££jj peak at x|z < 50 
or 65 < Mjj < 120GeV/c2, while all other back- 
grounds do not peak in both Z boson masses si- 
multaneously. We use a pythia Monte Carlo 
model [21] with the CDF detector simulation to 
estimate the small contribution from resonant di- 
boson processes and fit sideband data to collec- 
tively estimate all backgrounds that do not con- 
tain two bosons, collectively referred to as non- 
resonant background. 

We estimate the ££££ background by extrapo- 
lating the yield in the 185 < Mx < 300 GeV/c^ 
region to the signal region {Mx > 300GeV/c^ 
and Xzz ^ 5*-*) using a shape determined from 
a sample enhanced in non-resonant background. 
In order to construct samples enriched in this 
background, four-lepton candidates are selected 
in which some of the reconstructed leptons are 
"anti-selected" to fail one or more lepton identi- 
fication criteria. Anti-selected electrons must fail 
the HAD/EM energy selection and anti-selected 
muons must fail the minimum-ionizing energy se- 
lection. To further increase available statistics, 
the isolation requirement is removed for both cat- 
egories. Events reconstructed with the standard 
CDF processing that contain at least one trigger 
lepton and one anti-selected lepton are included 
in the reprocessing discussed in Section III. The 
Xzz ^'^x distributions for the resultant sam- 
ples are shown in Figs. 7 and 8. 

Figs. 5 and 6 show the invariant mass distri- 
butions of trigger lepton plus anti-selected lepton 
pairings for electrons and muons, respectively. 
The absence of an appreciable peak indicates 
few resonant Z boson events survive the anti- 



FIG. 5: Invariant mass distribution for pairs of can- 
didates consisting of a trigger electron and an anti- 
selected electron. 
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FIG. 6: Invariant mass distribution for pairs of can- 
didates consisting of a trigger muon and an anti- 
selected muon. 



selection. 

The two samples of four-body candidates 
that consist of a trigger lepton and either 
two- or three-anti-selected leptons with Mx > 
185GeV/c^ and Xzz < 500 are then fit simulta- 
neously to the empirical form 

to determine the falling shape of the Muu distri- 
bution (the power law parameter 7) and the rela- 
tionship of the number of events in the Xzz ^ 
ZZ window to the number in the off-mass side- 
bands (the exponential decay parameter r). As 
background composition and fake rate kinematic 
dependence varies with trigger and lepton type, 
we fit these sidebands separately for the eeee, 
eefj,fi, jijiee, and ^l^Ji^l^ji background shapes. Figs. 
9 through 12 show one-dimensional projections 
of the fit result for each channel against the fit- 
ted two- and three-anti-selected lepton x^z ^'^'^ 
Muu data as well as the one anti-selected sam- 
ple, which is not used in the fit. Table V lists 
the fit parameters obtained with their statistical 
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FIG. 7: Xzz vs. Afx distributions for the four-electron and electron-triggered two-electron two-muon sideband 
samples with 1, 2, and 3 anti-selected leptons. 



TABLE V: Four-lepton background fit results. 



Channel 


1 


r 


eeee 


-4.39 ± 0.09 


-0.0184 ± 0.0005 


ee/i/i 


-5.4 ± 0.2 


-0.0161 ± 0.0005 


/i/j,ee 


-5.3 ± 0.3 


-0.020 ± 0.002 




-6.5 ± 0.6 


-0.030 ± 0.003 



uncertainty. 

The background shapes obtained from these 
fits are normalized so that the sums of the in- 
tegrals for 185 < Mx < 300GeV/c2 and the 



simulation-derived diboson predictions match 
the number of events observed with 185 < Mx < 
300 GeV/c^ in the four-lepton samples. The 
shapes are then extrapolated into the low Xzz' 
high Mx m( signal regions. The statistical un- 
certainty on the normalization is the dominant 
source of uncertainty for the four-lepton non- 
resonant background prediction. 

As one test of the independence of the 
non-resonant predictions to the number of 
selected/anti-selected leptons, Tables VI and VII 
show the parameters and yield predictions ob- 
tained by fitting 7 for each sample independently 
of the others. All of the yield predictions for a 
given signal mass and decay channel are consis- 
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Xzz vs. Mx distributions for the muon-triggered two-electron two-muon and four-muon sideband 
samples with 1, 2, and 3 anti-selected leptons. 



tent with each other within the statistical uncer- 
tainty. 

The sideband data fit for the lijj non-resonant 
background estimates consist of events contain- 
ing a dilepton pair with x| < 25 and a dijet 
pair with either 40 < Mj, < 65 GeV/c^ or 120 < 
Mjj < 200GeV/c2. The Mjj spectrum near the 
Z boson pole mass is exponentially falling be- 
fore imposing any requirement on Mx but lin- 
ear for events with Mx > 300 GeV/ , where 
the effect of the four-body mass cut is to sculpt 
a peak in the dijet mass at Mjj ^ 200 GeV/c^ 
(see Fig. 13). We linearly interpolate the back- 
ground expectation for 65 < Mjj < 120 GeV/ 
from the lower and higher Mjj sideband data. 



To avoid underestimating the background at very 
high Mx where these sidebands are empty, we 
model the population of cither sideband vs Mx 
with an exponential fit to the available data to 
obtain the numbers used in the interpolation. 
Fig. 14 shows the numbers of events in the two 
dijet mass sidebands as a function of the require- 
ment on minimum four-body mass and the ex- 
ponential fits. Exponential functions model the 
data well. 

As one unbiased test of the prediction of the 
dijet mass spectrum, we repeat the selection and 
fit procedure on samples consisting of events con- 
taining a trigger Icpton plus an anti-selected lep- 
ton and at least two jets. Comparison of the fit 
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FIG. 9: Meeee and Xzz for ths one-, two-, and three- 
anti-selected (A-S) four-electron samples, and the re- 
sults of the simultaneous non-resonant background 
shape fit to the two- and three-anti-selected electron 
samples. 



FIG. 11: M„ 



and Xzz for the one-, two- 



and three-anti-selected (A-S) lepton samples for the 
muon-triggered two-electron two-muon channel, and 
the results of the simultaneous non-resonant back- 
ground shape fit to the two- and three-anti-selected 
lepton samples. 
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FIG. 10: Meenn and Xzz for the one-, two-, 
and three-anti-selected (A-S) lepton samples for the 
electron-triggered two-electron two-muon channel, 
and the results of the simultaneous non-resonant 
background shape fit to the two- and three-anti- 
selected lepton samples. 
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FIG. 12: M^^^p and Xzz for the one-, two-, and 
three-anti-selected (A-S) four-muon samples, and the 
results of the simultaneous non-resonant background 
shape fit to the two- and three-anti-selected muon 
samples. 



predictions against these 65 < Mjj < 120 GeV/c^ 
data, which are depleted of signal, show the 
method performs well (Tables VIII and IX) . The 
disagreement in the lowest M^ejj bin for this and 
other control samples is a result of a slight de- 
viation from an exponential distribution. This 
residual variation is taken as an extra systematic 
uncertainty in the lowest mass bin. 

We determine the backgrounds resonant in 
both Z boson masses with Monte Carlo mod- 
els normalized to the cross-sections predicted by 
MCFM [22]. In the Ujj channels, our dijet mass 
resolution and analysis selection does not distin- 
guish between the Z — >• jj££ and ZZ — >• jj£i, 
and so the background from both processes is 



present. 

Tables X through XV show the total prediction 
for each analysis channel. At each signal mass, 
the predictions are integrated over the four-body 
mass range listed in Table XVI. The uncertainties 
listed for the diboson predictions consist of the 
error on the mcfm cross-section 7%), the un- 
certainty on the luminosity (6%), and the statis- 
tical uncertainty due to finite Monte Carlo statis- 
tics, which is the dominant uncertainty on the di- 
boson prediction at high four-body mass, though 
a negligble component of the total background 
uncertainty. The uncertainties listed for the non- 
resonant backgrounds consist of the statistical 
uncertainty from the shape parameters and the 
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TABLE VI: Comparison of non-resonant background 
mass shape parameter 7 fitted independently in the 
individual sideband samples with the simultaneous fit 
to the two- and three-anti-selected Icirton samples. 



Channel 


eeee 








Simultaneous 


-4.39 


± 


0.09 


-5.42 


± 


0.15 


3 anti-leptons 


-4.22 


± 


0.11 


-5.19 


± 


0.17 


2 anti-lcptons 


-4.50 


± 


0.15 


-4.70 


± 


0.21 


1 anti-lepton 


-4.21 


± 


0.30 


-3.63 


± 


0.32 


anti-leptons 


-5.04 


± 


0.94 


-3.54 


± 


0.78 


Channel 










Simultaneous 


-5.25 


± 


0.34 


-6.51 


± 


0.61 


3 anti-leptons 


-4.98 


± 


0.50 


-6.3 


± 


1.1 


2 anti-leptons 


-4.96 


± 


0.41 


-6.60 


± 


0.73 


1 anti-lepton 


-5.25 


± 


0.73 


-5.33 


± 


0.92 


anti-leptons 


-5.7 


± 


1.6 


-4.4 


± 


1.3 



TABLE VII: Non-resonant background predictions 
for a characteristic example (the 410 590 GcV/c^ 
four-body mass range appropriate for a 500 GeV/c^ 
signal) from fits of the individual sideband samples, 
compared with the prediction from the simultaneous 
fit to the two- and three-anti-selected lepton samples. 
('liaiiiu>l (('((■ 



(Cflfl 



Simultaneous 
3 anti-leptons 
2 anti-leptons 

1 anti-lepton 
anti-leptons 



0.64 ± 0.29 
0.65 ± 0.31 
0.54 ± 0.26 
0.58 ± 0.30 
0.51 ± 0.35 



0.128± 0.064 
0.093± 0.048 
0.166± 0.091 
0.40 ± 0.24 
0.82 ± 0.59 



Channel 



Simultaneous 
3 anti-lcptons 
2 anti-leptons 
1 anti-lepton 
anti-leptons 



O.lSOtb 0.077 0.063± 0.040 

0.094t 0.060 0.058± 0.044 

0.127± 0.078 0.064± 0.043 

0.095± 0.065 0.17 ± 0.13 

0.14 ± 0.12 4.2 ± 3.8 



normalization uncertainty due to the small num- 
ber of events in the 185 < Mx < 300 GeV/c^ 
four-lepton control regions. The non- resonant 
and diboson background systematic imccrtaintics 
are negligible compared to the statistical uncer- 
tainty on the total background. 



VI. RESULTS 

We optimized all selections and estimated all 
backgrounds before examining the data with 
Mini > 300 GeV/ (? or with a dilepton pair hav- 



TABLE VIII: Comparison of dijet mass fit predic- 
tions with data for the trigger electron plus anti- 
selected electron and two-jet channel sample. Shown 
are the uncertainties on the mean prediction due to 

fit statistics. 



Meej-j (GeV/c^ 

350-450 
400-600 
500-700 
600-800 
650-950 
750-1050 
800-1200 



Prediction Observed 



1105 
488 
95.6 
18.8 
8.31 
1.64 ± 



± 29 
± 21 
± 7.2 
± 2.0 
± 1.1 
0.26 
0.73 ± 0.13 



941 
540 

115 
18 
8 
3 
2 



TABLE IX: Comparison of dijct mass fit predictions 
with data for the trigger muon plus anti-selected 
muon and two-jet channel sample. Shown are the 
uncertainties on the mean prediction due to fit statis- 
tics. 

M^njj (GeV/c^) Prediction Observed 



350-450 


23.7 


± 


4.0 


22 


400-600 


9.4 


± 


2.2 


9 


500-700 


1.48 


± 


0.59 


2 


600 800 


0.23 


± 


0.15 


1 


650-950 


0.093 


± 


0.072 


1 


750-1050 


0.015 


± 


0.017 





800-1200 


0.0059 ± 


0.0084 






TABLE X: Total eeee backgrounds with zz < 50 
for each signal mass Mx- The uncertainty includes 
the ZZ diboson production cross-section uncertainty, 
the luminosity uncertainty, the statistical uncertainty 
from the simulation, and the statistical uncertainty 
from the non-resonant background fit. 



Mx 
(GeV/c2) 


SM ZZ 


Non-resonant 


400 


0.22 ± 0.02 


1.31 ± 0.44 


500 


0.086 ± 0.009 


0.64 ± 0.29 


600 


0.045 ± 0.005 


0.44 ± 0.18 


700 


0.020 ± 0.003 


0.28 ± 0.20 


800 


0.007 ± 0.001 


0.15 ± 0.12 


900 


0.005 ± 0.001 


0.14 ± 0.11 


1000 


0.0032 ± 0.0001 


0.11 ± 0.10 
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FIG. 13: Dijet mass spectra for two-electron two-jet 
candidates with xL < 25 and three Meejj require- 
ments: no requirement, Meejj > 185 GeV/ c^, and 
300 < Meejj < 350 GeV/ c^, beyond which the shape 
of the dijet mass in the Z boson region is linear. 



TABLE XI: Total eefi^ backgrounds with x^zz < 50 
for each signal mass Mx- The uncertainty includes 
the ZZ diboson production cross-section uncertainty, 
the luminosity uncertainty, the statistical uncertainty 
from the simulation, and the statistical uncertainty 
from the non-resonant background fit. 



Mx 
(GeV/c2) 



SM ZZ 



Non-resonant 



400 
500 
600 
700 
800 
900 



0.19 

0.067 

0.035 

0.014 

0.004 

0.003 



± 0.02 
± 0.007 
± 0.004 
± 0.002 
± 0.001 
± 0.001 



0.33 ± 0.13 

0.128 ± 0.064 

0.075 ± 0.047 

0.041 ± 0.029 

0.019 ± 0.013 

0.016 ± 0.013 



1000 0.0013 ± 0.0006 0.012 ± 0.011 
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FIG. 14: Number of two-lepton two-jet events with 
40 < Mjj < 65 GeV/c^ and 120 < Mjj < 200 GeV/c^ 
as a function of the minimum Mujj mass require- 
ment, along with the exponential fit to each that is 
used to interpolate the background prediction for the 
65 < Mjj < 120 GeV/c^ region. 



TABLE XII: Total /ifiee backgrounds with x^zz < 
50 for each signal mass Mx- The uncertainty in- 
cludes the ZZ diboson production cross-section un- 
certainty, the luminosity uncertainty, the statistical 
uncertainty from the simulation, and the statistical 
uncertainty from the non-resonant background fit. 



ing Xz < 25 and a dijet pair with 65 < Mjj < 
120 GeV/ c^. Figs. 15 and 16 show the data 
in these regions and the combined resonant and 
non-resonant background predictions for all four- 
lepton channels and for both dijet channels. In 
all cases the data agree with the total back- 
ground prediction and provide no compelling ev- 
idence for resonant ZZ diboson production. The 
highest-mass ££M event (577 GeV/c^) consists of 
four muons. For this event, one Z boson candi- 



Mx 
(GcV/c^) 



SM ZZ 



Non-resonant 



400 0.077 ± 0.008 

500 0.027 ± 0.003 

600 0.014 ± 0.002 

700 0.0065 ± 0.0010 

800 0.0018 ± 0.0007 

900 0.0014 ± 0.0006 

1000 0.0011 ± 0.0005 



0.32 ± 0.16 
0.130 ± 0.077 
0.078 ± 0.055 
0.044 ± 0.034 
0.021 ± 0.017 
0.018 ± 0.017 
0.014 ± 0.013 
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TABLE XIII: Total fifififi backgrounds with x^zz < 
50 for each signal mass Mx- The uncertainty in- 
cludes the ZZ diboson production cross-section un- 
certainty, the luminosity uncertainty, the statistical 
uncertainty from the simulation, and the statistical 
uncertainty from the non-resonant background fit. 



2.5-2.9 fb ' 



M 



X 



(GeV/c2) 



SM ZZ 



Non-resonant 



400 
500 
600 
700 
800 
900 
1000 



0.090 

0.036 

0.018 

0.0082 

0.0018 

0.00011 

0.0009 



± 0.010 0.21 

± 0.005 0.063 

± 0.002 0.031 

± 0.0015 0.015 

± 0.0007 0.0056 

± 0.00005 0.0046 

± 0.0005 0.0031 



± 0.11 
± 0.040 
± 0.023 
± 0.013 
± 0.0049 
± 0.0042 
± 0.0030 



C = 2.5-2.9 fb-i 



m 



> IQ-l - 




400 500 600 700 800 900 1000 1100 1200 

M„„ (GeV/cO 

FIG. 15: Prediction and data for all four-lepton 
channels combined. The background prediction for 
each bin consists of the integral of the non-resonant 
background functions and diboson Monte Carlo de- 
termined in Section V. The background predictions 
agree with the data. 



date has a mass of 79±4.2 GeV/c^ and the other Z 
boson candidate has a mass of 400 ± 170 GeV/c^. 
The large mass uncertainty of the second Z bo- 
son comes from a large curvature uncertainty in 
the measurement of one px = 290 GeV/c muon 
with few COT hits. The highest-mass iijj event 
(868 GeV/c2) has M^e = 96.5 ± 1.3GeV/c2 and 
Mjj = 77.8 ± 6.5 GeV/c2. 

Absent evidence of a signal, we communicate 
our sensitivity to X ZZ processes by set- 
ting limits with an acceptance from a widely- 
available HERWIG Monte Carlo process [23], the 
spin-2 Kaluza-Klein graviton. The total accep- 
tance times efficiency for this process varies be- 
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FIG. 16: Prediction and data for both Ujj chan- 
nels combined. The background prediction for each 
bin consists of the integral of the non-resonant back- 
ground functions and diboson Monte Carlo deter- 
mined in Section V. The background predictions 
agree with the data. 



tween roughly 40-50% for a given four-lepton 
channel and between 20-40% for a given dijet 
channel (Fig. 17). At higher masses, the fraction 
of Z boson decays increases in which the angu- 
lar separation between products is too small for 
the calorimeter to resolve. This lowers the accep- 
tances for the dijet modes and, to a lesser extent, 
the electron modes. 

The combined effect of the lepton reconstruc- 
tion and identification improvements on graviton 
signal is demonstrated in Figs. 18 and 19 for the 
four-electron and four-muon channels. We com- 
pute data yields and estimates of signal and back- 
ground by integrating the Monte Carlo predic- 
tions and fitted non-resonant background shapes 
over a set of overlapping, variable-width bins for 
signal masses from 400 GeV/c^ to 1 TeV/c^ (Table 
XVI) . Each signal bin width is chosen to be large 
enough to fully contain the four-body mass distri- 
bution expected for an intrinsically narrow signal 
and the broadening from systematic effects. Ta- 
ble XVII shows the total background prediction 
and observed data yields in each of these bins. 

We calculate 95%-credibility upper limits as a 
function of signal mass using a six-channel prod- 
uct of Bayesian likelihoods and a uniform prior 
for the (non-negative) X ZZ cross-section. 
We use marginalized truncated-Gaussian nui- 
sance parameters for the luminosity, background 
predictions, and signal efficiencies, and we ac- 
count for systematic uncertainties correlated 
amongst the six channels when appropriate. As 
discussed earlier, we assign a 20% uncorrelated 
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TABLE XIV: Total eejj backgrounds with 65 < Mjj < 120GeV/c^ for each signal mass Mx- The uncertainty 
includes diboson cross-section uncertainties, the uncertainty on the luminosity, the statistical uncertainty from 
the simulation, and the uncertainties from the non-resonant background fits. 



(GeV/c2) 



SM ZZ SM W^Z Non-resonant 



400 5.72 ± 0.97 9.4 ± 1.1 483 ± 18 

500 2.43 ± 0.58 3.25 ± 0.46 128.0 ± 8.2 

600 0.99 ± 0.36 1.10 ± 0.22 47.4 ± 4.1 

700 0.19 ± 0.18 0.60 ± 0.16 14.9 ± 1.7 

800 0+° " 0.158 ± 0.083 2.86 ± 0.46 

900 0+0 " 0.095 ± 0.067 1.75 ± 0.31 

1000 0+° " 0+° °*^^ 0.77 ± 0.16 



TABLE XV: Total nfJ.jj backgrounds with 65 < Mjj < 120GeV/c^ for each signal mass Mx- The uncertainty 
includes diboson cross-section uncertainties, the uncertainty on the luminosity, the statistical uncertainty from 
the simulation, and the uncertainties from the non-resonant background fits. 



Mx 
(GeV/c2) 


SM ZZ 


SM 


Non-resonant 


400 


2.90 ± 0.57 


6.04 ± 0.73 


162 ± 11 


500 


1.30 ± 0.38 


2.06 ± 0.32 


37.7 ± 4.4 


600 


0.57 ± 0.26 


0.73 ± 0.17 


12.6 ± 2.0 


700 


0.26 ± 0.19 


0.229 ± 0.93 


3.53 ± 0.72 


800 


0.09 ± 0.13 


0.023 ± 0.040 


0.57 ± 0.16 


900 


Q-(-0.10 


Q+0.032 


0.33 ± 0.10 


1000 


Q-(-0.10 


Q+0.032 


0.133 ± 0.045 
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FIG. 17: Products of acceptance times efficiency for 
each of the graviton analysis channels and their de- 
pendence on graviton mass. These do not include ZZ 
diboson branching ratios, which for each ZZ — > ££jj 
mode are approximately 40 times the branching ra- 
tios for ZZ — >■ eeee or ZZ — >■ ^/x/i/i. The ee/i/i and 
Hfiee acceptances have been summed. 
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FIG. 18: Four-electron yield comparison for a 500 
GeV/ graviton between (lower histogram) the CDF 
standard electron selection criteria and (upper his- 
togram) the criteria employed in the present analy- 
sis. 



uncertainty to the total acceptance x efficiency 
for each channel to account for the observed time- 
dependent variation in the DrcU-Yan fifi cross- 
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FIG. 19: Four-muon yield comparison for a 500 
GeV/c^ graviton between (lower histogram) a CDF 
standard muon selection criteria and (upper his- 
togram) the criteria employed in the present analysis. 



TABLE XVI: Signal binning used for limit-setting. 
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FIG. 20: 95%-credibility cross-section upper limit 
assuming acceptance for a massive graviton along 
with the limit and 68% variation expected for the 
background-only hypothesis. 



Signal Mass 
(GeV/c^) 



Bin Half 
Width (GeV/c2) 



400 ± 
500 ± 
600 ± 
700 ± 
800 ± 
900 ± 
1000 ± 



70 
90 
130 
160 
160 
230 
280 



ground estimates derived from sideband data 
samples and electroweak Monte Carlo models. 
To quantify our sensitivity, we set limits using 
the acceptance for a Randall- Sundrum graviton 
model that are 7-20 times stronger than the pre- 
viously published direct limits on resonant ZZ 
diboson production. 



section, conservatively covering the sum of in- 
dividual systematic uncertainties such as signal 
acceptance uncertainties in order to simplify the 
combination. Studies of the individual uncertain- 
ties indicate the largest contribution after the un- 
certainty due to the Z — >■ /i/i cross section vari- 
ation is the 5.9% uncertainty on the luminosity. 
In addition to the observed limit, we compute ex- 
pected limits from 10,000 pseudo-experiments at 
each candidate X mass. Fig. 20 shows the resul- 
tant limits along with the k/Mpi —0.1 Randall- 
Sundrum (RSI) graviton cross-section from her- 
WIG. The present search improves the 0(4 pb) 
limit of the earlier eeee search [9] by an order of 
magnitude. 



VII. CONCLUSIONS 

We have reported on an improved search for a 
massive resonance decaying to ZZ dibosons via 
the eeee, eefifj,, fJ-iJ-fJ-fJ-, eejj, and /i/ijj channels. 
We find that the four-body invariant mass spec- 
trum above 300 GeV/c^ is consistent with back- 
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TABLE XVII: Total background prediction and observed data yields for each of the limit-setting bins in Table 
XVI. Successive bins are partially correlated. The uncertainty (quoted as the least two significant figures in 
parentheses) is the systematic uncertainty on the mean background prediction and does not include statistical 
fluctuation about the mean. 



Channel 


400 


Mx (GeV/c2) 
500 600 700 800 


900 


1000 


eeee 


Expected 1.53(44) 
Observed 


0.73(29) 0.49(18) 0.30(20) 0.16(12) 



0.15(11) 



0.11(10) 



eefifj, 


Expected 0.52(13) 
Observed 2 


0.195(64) 0.110(47) 0.055(29) 0.023(13) 



0.019(13) 



0.013(11) 



/i/iee 


Expected 0.397(16) 0.157(77) 0.092(55) 0.050(34) 0.023(17) 
Observed 


0.019(17) 



0.015(13) 





Expected 0.30(11) 
Observed 1 


0.099(40) 0.049(23) 0.023(13) 0.0074(49) 0.0047(42) 0.0040(30) 
1110 


eejj 


Expected 498(18) 
Observed 456 


133.7(82) 49.5(41) 15.7(17) 3.02(48) 
142 69 28 7 


1.84(34) 
5 


0.77(21) 
2 




Expected 171(11) 
Observed 143 


41.1(44) 13.9(20) 4.02(75) 0.68(21) 
41 19 4 2 


0.33(14) 
2 


0.13(11) 
1 



APPENDIX 

The standard CDF reconstruction software 
uses two main approaches to reconstruct tracks. 
High quality central tracking (I77I < 1) starts in 
the COT and assembles picccwisc segments of 
up to 12 hits in each supcrlayer, fits them, and 
groups them into tracks to which any available 
silicon hits are then attached in an outside-in 
search. Afterward, "silicon standalone" track- 
ing starts with all possible combinations of three 
unused silicon hits, searches the remaining sil- 
icon layers, and projects successful tracks into 
the COT to attach any compatible hits in order 
to improve the track momentum resolution and 
lower the fake rate. 

The combination of these approaches results in 
low efficiency in the 1 < jr;! < 2 region. Tracks 
originating from z — Q with |r/| < 1.7 will leave 
traces of their passage in the lowest-radii super- 
layers of the COT. Though very efficient when 
full COT coverage is available, for I77I > 1 the 
central tracking algorithms lose efficiency nearly 
linearly with |ry| reaching zero efficiency at about 
|?7| = 1.6. The silicon fully covers \rj\ < 1.8 to 
compensate for the falling COT efficiency, but 
the existing silicon-driven tracking algorithms re- 
construct tracks with low efficiency and produce 
low-quality or spurious tracks with poor pointing 
resolution into the COT. Thus the COT informa- 
tion for forward tracks is rarely exploited. 



This analysis employs a thorough revision of 
the forward and central tracking algorithms in 
order to reconstruct tracks with better efficiency 
and resolution, including a new "Backward" al- 
gorithm that makes full use of the partial COT 
coverage. The Backward algorithm, illustrated 
in Fig. 21 for a simple case, starts by search- 
ing the COT for hits unused by the central COT 
algorithm and constructing segments in one of 
the inner axial superlayers consisting of no more 
than 12 hits. At this stage, the position mea- 
surements contain a drift sign ambiguity and im- 
portant drift time corrections, such as large time 
of flight and sense wire signal propagation times, 
are unknown and cannot be approximated by the 
constant corrections assumed for the central seg- 
ment pattern recognition. The Backward algo- 
rithm solves this problem with a variant of the 
central segment pattern recognition that resolves 
the drift sign ambiguity and drift time correc- 
tions during the search and is optimized for track- 
ing in the low radius, high hit density inner su- 
perlayers and near the COT endplates. Once 
unused COT hit segments are found consistent 
with a forward track, the algorithm then fits the 
segments with a beamline constraint to obtain 
five-parameter helices that intersect the z posi- 
tion of the highest sum z vertex identified 
using central algorithm tracks. In most cases, 
the fits do not conclusively identify stereo COT 
measurements in the innermost stereo superlayer, 
and so multiple helices are obtained correspond- 
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FIG. 21: A simple example of the Backward tracking 
algorithm in low luminosity data. In a stretched and 
rotated r — <j) view of the relevant section of the track- 
ing volume, COT pulses in a single 12-layer axial su- 
perlayer indicating two possible hit locations corre- 
sponding to —(j) and +<!} drift are processed to identify 
trajectory segments and fitted to obtain drift time 
corrections and an initial trajectory with large un- 
certainties. An iterative Kalman filter search through 
possible rj values for silicon charge clusters consistent 
with the initial trajectory produces a tree of track 
possibilities, from which the single best candidate, 
shown with a projection of the final Scr uncertainties, 
is chosen. Also shown is an independent measure- 
ment from the forward calorimeter shower maximum 
scintillator. 



ing to trajectories through each possible combi- 
nation of sihcon module and polar angle. Drift 
time corrections are recomputed for each case. 

After obtaining initial helix fits, including esti- 
mates of the helix uncertainties and correlations 
obtained from the small number of COT hits, the 
algorithm begins an outside-in silicon hit search 
that uses a Kalman filter to correct for energy 
loss and multiple scattering in the tracker mate- 
rial. After each search completes, quality criteria 
based on hit pattern, multiple usage, charge de- 
position, and module operational status are ap- 
plied to the successful hypotheses, and all surviv- 
ing hypotheses except those with the most hits 
are discarded. Finally, the COT is searched for 
any remaining information, including stereo mea- 
surements in inner superlayers. 

The Backward algorithm has been validated 
on a variety of samples, with emphasis on large 
samples of Z — > ee and Z — > fifj, simulation 
and data. Fig. 22 shows the improvement in 
Z — > fxfi yield in muon-triggered data involving 
higher-quality forward tracks with COT hits in a 
subset of the data, demonstrating the increase in 
muon acceptance due to the new software. The 
lower curve represents the dimuon mass spectrum 
for the combination of a trigger muon tracked 
with the central algorithm and a forward muon 
tracked in the combination of the COT and the 
silicon detectors with the silicon-driven algorithm 
in the standard software. The upper curve shows 
the same spectrum in the new software, where 
the Backward algorithm has largely superseded 
the other silicon-driven algorithm. With a mod- 
est increase in background, the peak yield has 
improved by about 260%, corresponding to an 
approximately 10% increase in the total Z iifi 
yield over the entire detector. The distributions 
of all forward muon identification variables are 
qualitatively the same as those of muons found 
with the central COT-driven algorithm, indicat- 
ing that we have selected a sample of forward 
muons with purity comparable to the central 
muons. 
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